Emergence of a common energy scale close to the orbital-selective Mott transition 
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We calculate the spectra and spin susceptibilities of a Hubbard model with two bands having 
different bandwidths but the same on-site interaction, with parameters close to the orbital-selective 
Mott transition, using dynamical mean-field theory. If the Hund's rule coupling is sufficiently strong, 
one common energy scale emerges which characterizes both the location of kinks in the self-energy 
and extrema of the diagonal spin susceptibilities. A physical explanation of this energy scale is 
derived from a Kondo-type model. We infer that for multi-band systems local spin dynamics rather 
than spectral functions determine the location of kinks in the effective band structure. 

PACS numbers: 71.27.+a, 71.30.+h 



The interactions in correlated metals lead to the emer- 
gence of characteristic energy scales. Close to the Fermi 
energy Landau Fermi-liquid theory pQ applies and the 
effective electronic dispersion E^ is renormalized, but 
remains linear as in the noninteracting case. The lin- 
ear dependence terminates at an excitation energy which 
cannot be calculated within Fermi-liquid theory itself. 
With increasing strength of the electron-electron interac- 
tion this Fermi-liquid coherence scale decreases and ul- 
timately vanishes at the Mott transition from a metal 
to an insulator. At the same time charge excitations 
are shifted to higher energies of the order of the in- 
teraction energy and are thus suppressed. For the 
single-band Hubbard model with on-site Coulomb inter- 
action U it was shown [2 in dynamical mean-field theory 
(DMFT) [3] that the Fermi-liquid regime terminates at 
an energy scale ws at which the real part of the self- 
energy, Re[£(cj)], and hence the effective dispersion E^, 
has a rather sudden change in slope [I]. This 'kink' does 
not require any coupling to external bosonic degrees of 
freedom but is due to the correlated behavior of interact- 
ing electrons. For the single-band model the Fermi-liquid 
scale loy. can be derived from the low-energy properties of 
the local spectral function [2]. Moreover, it was demon- 
strated |3] that the energy scale ws is linked to the 
characteristic energy scale w sp of spin fluctuations. Kinks 
in the electronic dispersion were studied theoretically in 
a variety of contexts |7Hl4] . 

In this Letter we explore the origin and characteristic 
energy scale of kinks in the effective electronic disper- 
sion in a more general context. Employing the DMFT 
we study a two-band Hubbard model with two differ- 
ent bandwidths, the same on-site repulsion U for both 
bands, and an interorbital repulsion U\ and ferromag- 
netic Hund's rule spin exchange J between the bands. 
Thereby it is possible to capture orbital effects in corre- 
lated materials that do not exist in single-band models. 
Indeed, different kinks in the dispersion depending on the 
orbital character are observed for S^RuO^ both experi- 
mentally and theoretically [T51 ITB] . 



We study the model Hamiltonian 

ijma 

flint = U ^ ni m ^n im ± + / ](Ul ~ Saa'J) 7W™i2o-' 
im iaa' 

+ ^ ^imai^ima^ima + ^\ma^ima)^ima ' 

with spin index a = f , 1 and orbital index m = 1, 2. Here 
a bar over an index denotes the opposite spin or orbital. 
The two bands do not hybridize but are coupled by the 
interorbital interactions U\ and J. We consider U, U\ 
and J as independent parameters which can take arbi- 
trary values, but U\ = U — 2J [T7] for d electrons. As 
in the single-band case the correlation strengths of the 
two bands may be roughly parametrized by the ratios 
U/Wi and U/W2, which are assumed to be unequal. Due 
to this difference in the relative interaction strengths an 
orbital-selective Mott transition (OSMT) occurs upon in- 
crease of U |18ff3"3"] . We assume semi-elliptic densities of 
states, p m (e) = (8/n)^{W m /2) 2 - e 2 /W^, with band- 
widths Wi < W2. Since the hopping amplitudes Uj. m 
are diagonal in the band index, so are the single-particle 
Green functions and self-energies. Off-diagonal contribu- 
tions only occur in two-particle and higher-order correla- 
tion functions, e.g., spin and charge susceptibilities. 

In DMFT the model (JT|) is mapped onto the following 
two-impurity Anderson model (TIAM) |28[ [33] 

#TIAM = ^2 £kmC krna C k,na + E ^^ma 
kma ma 

+ E {y^lmJma + ^.C.) + , (2) 

kma 

where the local interaction H\°£ has the same form 
as H int , but without the index i. The DMFT self- 
consistency conditions demand that the band energies 
ekm and hybridizations Vkm are determined such that 
Green functions and self-energies of (p| equal the corre- 



2 



1.2 



0.x 



0.4 



.(a) 




A! 


Wi — 


0.4 


W 2 


= 1.4Wi 


A 2 


Wi — - 


0.0 


J 


= 


U 


= 1.4Wi 




-Ui 


= u / 2 I 






-0.4 










1.0 








V 


0.0 










-1.0 





\ ReZi/Wi — ■ 




X Re£ 2 /Wi — 














. I Y S i P -W! — 




J^/j 6 \2 ■ Wi — - 



1,-, 







-0.125 





j/Wi 



0.125 



FIG. 1. For J — and U\ ^= the difFerence in the correlation 
strength of the two bands is observed not only in the shape of 
the spectral function A m (oj) (a), but also in the corresponding 
band-resolved Kondo temperatures and widths of the Fermi- 
liquid regime (b, b'). 



sponding local lattice quantities. 

Pm(e) 



G m (w) = de 



(3a) 



'u + iO - e m - S m (w) - e ' 
£m<rM = cj + iO - e m - 1 - A m (w) , (3b) 

which take the same form as for two decoupled one- 
band models due to the absence of interorbital hopping. 
Here the hybridization function is defined as A m (uj) = 
J2k l^fcm| 2 /( w + *0 — e km ). We solve the impurity model 
using the Numerical Renormalization Group (NRG) . The 
complete Fock space NRG [33] proceeds as in the single- 
band Hubbard model |35| . but with a local dimension 
of 16 for the impurity and the chain sites |19[ [35] . Un- 
less noted otherwise the NRG discretization parameter 
is A = 2.5, and we keep on the order of 10 5 states, in- 
cluding multiplicities of irreducible subspaces, in each 
NRG iteration. Although we focus on the low-energy 
region we employ Oliveira's [36 z-trick (N z — 4) to im- 
prove the spectra at higher energies. To obtain the self- 
energy with high quality the correlation function F m (uj) 
— ((Mm <r! ^intli ^m,o-))w i s calculated since it is numeri- 
cally better conditioned than the Dyson equation [37 . 

We compute the spectra, self-energies, and spin sus- 
ceptibilities for the two bands in the metallic phase close 
to the OSMT, and monitor the behavior as a function 
of the Hund's rule coupling J. The different correla- 
tion strengths of the orbitals lead to different behavior of 
the spectral functions A m (uj) = — Im[G m (w)]/7r and self- 
energies E m (w) which are shown in Figs. [T] and [2] The 
overall behavior corresponds to that of two Fermi liq- 
uids with different mass renormalizations, reminiscent of 
two uncoupled one-band Hubbard models with different 
local interactions. In all cases the single-particle spec- 
tra Ai(uj) and ^(w) differ significantly, especially close 
to the OSMT when the spectral function of the narrow 
band has a very sharp central peak. Because the spec- 
trum of each band depends mostly on the correlation 
strength U/W m but not much on Hund's rule exchange 
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FIG. 2. In the metallic phase close to the OSMT a striking 
proportionality, Re[Ei] oc Rep2] and Xi P ^ X2 (panels b, 
b', d, d') is seen to emerge for J > 0, which is in marked 
contrast to the decoupled behavior in Fig. [T] Very close to 
the OSMT (d, d') even a weak exchange coupling J leads to 
a common low-energy scale. Note that U\ has little influence 
on the qualitative low-energy behavior. 



J, the latter was previously characterized 'band 
decoupler' [26J, at least for weak interorbital hopping. 
Quantum-Monte Carlo results [22 123 [23 [551 \M suggest 
that at T = the low-frequency behavior of Re[£ m ] and 
Xm is then also different. Fig. [T] shows that this is indeed 
the case — but only for J — and U\ =^ 0. Indeed, a 
finite Hund's rule coupling leads to a fundamentally dif- 
ferent low-energy behavior of Re[£ m ] and Xm- Namely, 
as the system approaches the OSMT we find that at 
low energies these quantities become proportional, i.e., 
Re[Ei(w)] cx Re[E 2 (w)] and x ?{u) cx x?^) EE]- As 
illustrated in Figs. [2})', [2ji, [2}T this striking result 
cannot be inferred from the spectral functions A m (uj) 
since the shape and the characteristic energy scales of 
the latter differ considerably and thus suggest a decou- 
pled behavior. The characteristic energy scale of the spin 
fluctuations, i.e., the locations of the extrema in the 
spin susceptibilities Xmi 1 ^)- allow one to define a Kondo 
temperature for each band. The proportionalities dis- 
cussed above imply that the system has identical Kondo 
temperatures (w^ p = to^), Fermi-liquid energy scales and 
kinks (uif = w|"), irrespective of the different correla- 
tion strengths of the bands. Furthermore, the self-energy 
kinks and the strongest spin fluctuations occur at the 
same energy in each band, ~ to^ (as observed also 
in the single-band case [5]), which means that for the 
two-band model (JTJ a single common low-energy scale 
emerges for kinks and spin fluctuations in both bands. 
For our numerical data we define the kink scale w£ as 
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FIG. 3. Intensity plot of the total momentum-resolved spec- 
tral function Afe(w) deep inside the central peaks as a function 
of efc and lo for the same parameters as in Fig. |2^i,b. The ef- 
fective dispersion 75*., m is given by the local maxima in Ak(ui) 
(blue line for the wide band, lighter red line for the narrow 
band). It is linear (dashed lines) in the Fermi- liquid regime 
close to the Fermi surface and has kinks at the same energy 
~ ±0.0028W / i for both bands (solid horizontal lines). 

the energy for which the extrapolated linear dispersion 
near the Fermi energy deviates from Re[E m (w)] by 20%, 
which agrees well with the perceived location of the kinks 
in Figs. [2|),d. The corresponding momentum-resolved 
spectral function Ak{oj) and effective dispersions Ek m 
are shown Fig. [3j We observe that although the slope of 
the Fermi-liquid dispersion is very different for the two 
bands, the linear regimes terminate at the same energy 
scale, which however slightly deviates from due to 
band structure effects. 

By comparing the results for J = in Fig. [I] with 
,7^0 in Figs. [2Jd-c it is clear that the interorbital repul- 
sion Ui is not responsible for the common energy scale. 
This effect only appears in the presence of the Hund's 
rule coupling J, whereas U\ merely leads to quantitative 
modifications. We will thus restrict ourselves to U\ = 
in the following. Starting from J — we study the con- 
tinuous evolution of the two initially uncoupled Hubbard 
models into the 'locked' regime. To this end we obtain 
Kondo scales to^(J) and self-energy kinks w^(J) for the 
two orbitals for two different values of U /W\ (Figs.[4ji-b). 
As J is increased both orbital-resolved energy scales ap- 
proach each other and finally merge into a single scale, as 
seen in Figs.[2j>d'. Comparing Figs.[4jt, and b, we observe 
that this common low-energy scale appears at a thresh- 
old value which decreases for increasing U. We also no- 
tice the very close correspondence between kink energies 
and Kondo-temperatures, especially in the more strongly 
correlated case (Fig.jij}). This observation can be under- 
stood in terms of local Fermi-liquid theory Since 
the binding energy of the Kondo singlet is approximately 
given by the Kondo temperature, the linear regime must 
terminate at ui ~ w sp ; see Refs. [6] for a discussion 
of the single-band case. As expected the kink energy 
scale derived in Ref. [2J, Lo* n — 0.2Z m W m , applies only 
to the narrow band with its well-developed three-peak 
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FIG. 4. (a,b) Orbitally resolved kink energies ujf 2 (J) and 
Kondo temperatures ^ P 2 (J) as calculated within DMFT for 
different Hubbard interactions U. Sufficiently large U leads 
to [^^(J) = cj^(J) and sufficiently large J leads to cof(J) = 
oj^(J) and ul P (J) = u^'(J). The dashed lines in (b) mark 
the single-band estimate for [2] which applies only to the 
narrow band (see text), (c) A simplified Anderson impurity 
model and, (d), a related two-impurity Kondo model, both 
with a behavior similar to the DMFT solution. Note that 
the system in (b) enters the OSMT phase at approximately 
J — 0.2(7. For these very low-energy features 0(5500) states 
were kept and N z — 2. 



spectral function and not to the wide band (cf. solid and 
dashed lines, respectively, in Fig. [4J3, corresponding to 

o.s< m)- 

In the vicinity of the OSMT the characteristic ener- 
gies u^(J) and oj^(J) represent equivalent energy scales 
and hence contain the same physical information. In or- 
der to explain the J dependence it suffices to discuss one 
of them, and we will focus on uj^(J) in the following. 
To explain the locking of the low-energy scales for the 
two bands we proceed in two steps; see Figs. [2J: and[3Ji. 
First we establish that the locking is an intrinsic prop- 
erty of the underlying TIAM Hamiltonian and is only 
quantitatively modified by the DMFT self-consistency 
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equations Q3J) - Then we compare with the results for a 
Kondo-type model that allows us to identify the compet- 
ing couplings and elementary excitations. For the first 
step we solve the impurity model ^ with different but 
constant hybridization functions (A2(w) = 1.4Ai(cj) = 
const) and extract the Kondo temperatures ut^(J) from 
the maxima of the spin susceptibilities. The result is 
depicted in Fig. [4J3 for two values of U, showing very 
good qualitative agreement with the DMFT results in 
Figs. [4^, [4Jd. In particular, the common low-energy scale 
emerges at a value of J which decreases with increasing 
U in a similar fashion. We conclude that the DMFT self- 
consistency induces only minor modifications as long as 
the system remains in the metallic phase. 

In a second step we focus on the low-energy spin dy- 
namics close to the OSMT. In this regime charge exci- 
tations are strongly suppressed. We thus consider the 
Kondo limit of ([2| which captures the low-energy spin 
dynamics of the TIAM in Fig. [4};. In this limit the 
Hamiltonian (I2J) reduces to a two-impurity Kondo model 

(2ikm) unrig, 

H2IKM = e kmcl ma C km(J 
kma 

J5i-S 2 . (4) 

m 

Here J > is the Hund's exchange interaction of 
while the antiferromagnetic couplings J m stem from su- 
perexchange processes and decrease with increasing U. 
We take J2P2(0) = 1.4 Jipi(O) to obtain different Kondo 
temperatures for J = 0, i.e., w^ p (0) ^ u^(0). The J de- 
pendence of uj^(J) is shown in Fig. The qualitative 
agreement among the results obtained for all three mod- 
els (DMFT, TIAM with constant hybridization, 2IKM) 
confirms that Q already describes the essential processes 
that lead to the emergence of the joint low-energy scale. 
In the 2IKM the spins will align for low-excitation ener- 
gies and form a composite spin-1 object |41) . This hap- 
pens roughly when the energy gain ui s 1 p ( J) + lu^(J) due 
to Kondo screening of the two impurities is overcome 
by J/4, the approximate energy gain due to the ferro- 
magnetic exchange. Hence the locking of the low-energy 
scales sets in at about J « w^p(J)/8, as seen in Fig. [4^ 
(J w 0.2Wi) and|4jD (J w O.lWi). We conclude that the 
low-energy spin dynamics of the two impurity spins (and 
thus the two bands of the corresponding lattice model 
in DMFT) exhibit joint fluctuations and thus have equal 
Kondo scales if J dominates the individual Kondo scales 
wfP and are essentially independent otherwise. Regard- 
ing the influence of U, we note that the antiferromagnetic 
couplings between the spins and the baths decrease with 
increasing interaction, i.e., more correlated systems ex- 
hibit stronger locking of their spins and their low-energy 
scales. Thus J couples the low-energy scales more effec- 
tively for stronger correlations, as seen in Fig. [4jt-c. 
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FIG. 5. The correlation functions xtpi^) an d Xsp( w ) indicate 
the creation of a composite spin-1 object for large J (green 
curves in (a) and in the inset have about the same magnitude; 
also note the different scales in (a) and the inset). Accord- 
ingly, the partial spectral weight fraction Xpart./Xp ar t. m ( b ) 
grows strongly with J and is only quantitatively affected by 
the upper limit u. Here uio = max[cj° p (J),uj2 P (J)]. 

To explicitly verify the physical picture described 
above we investigate the correlation functions 

X t p (uj) = ((S 1 ±S 2 \S 1 ±S 2 )) u (5) 

in DMFT. They describe the dynamics of the composite 
spin-1 object and the 'residual' singlet x~(w), re- 

spectively, and are plotted in Fig. [5ji. As expected, for 
both susceptibilities the positions of their maxima de- 
creases with increasing interaction. However, only the 
composite spin-1 object shows a resonance that also in- 
creases in amplitude. In Fig. |5]d we plot the integrated 

weight Xpart — Jo X ± (<*0^ w as a function J for several 
values of the upper limit w. In particular, for larger val- 
ues of J the residual spin contributes only little to the 
total low-energy spin response xT + X2 P °f ^ ne system, 
which is thus well described by x + ( UJ )- This establishes 
the formation of the composite spin-1 object as the physi- 
cal origin for the emergence of the common energy scale. 
Namely, as the OSMT is approached by increasing U, 
the antiferromagnetic superexchange of the narrow and 
wide band both decrease until J becomes the dominat- 
ing scale, at least for the spins in the narrow band. The 
spins align and create a composite spin-1 and exhibit 
joint low-energy dynamics, leading to the proportionali- 
ties of Re[£ m (w)] and Xmi 00 ) f° r * ne two bands. 

In summary, we explored the physical mechanism for 
kinks responsible for the appearance of kinks in the 
self-energy by studying a two-band model with differ- 
ent bandwidths but the same local charge interactions as 
well as Hund's rule spin exchange. We find that the phys- 
ical picture developed previously for single-band systems 
close to the Mott transition is significantly modified for 
strong Hund's rule coupling, due to the formation of a 
local spin-1 object. As a consequence, a common low- 
energy scale emerges for the two bands for the kinks in 
the self-energies and the maxima in the spin susceptibil- 
ities. 
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